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INTRODUCTION 


1.0 

1.1  BACKGROUND 

The  subject  of  extending  survey  control  into  unchartered  areas  has 
been  of  interest  for  many  years.  The  AN/USQ  28  system  was  one  of  the  more 
well  known  systems.^ ) As  progress  in  mapping  continues  the  efficiency  of 
such  systems  will  continue  to  improve  but  there  is  continuing  need  to 
explore  unique  ways  in  which  to  implement  such  systems. 

Analytic  aerotriangulation  has  advanced  very  rapidly  in  the  past 
decade  and  under  ideal  conditions  it  can  now  compete  with  geodetic  surveys 
in  the  determination  of  survey  coordinates  for  objects  that  are  imaged  on 
the  photograph  of  interest.  A certain  number  of  control  points  are  needed 
in  order  to  perform  the  aerotriangulation.  In  most  cases  these  objects 
are  stationary  points  located  at  ground  level.  The  selection  of  control 
points  and  the  determination  of  these  positions  can  take  place  in  several 
ways.  Basically  they  are  determined  either  before  the  photography  or  after 
the  photography.  Each  procedure  has  its  advantages.  In  either  case  the 
flexibility  of  control  location  is  limited. 

The  concept  of  a moving  airborne  control  target  offers  a combination 
of  the  advantages  of  each  control  procedure  mentioned  above.  The  airborne 
control  target  can  be  easily  seen,  it  can  be  positioned  as  necessary, 
conceptually  we  can  have  as  many  airborne  control  images  as  we  please,  they 
will  not  be  destroyed  or  moved,  and  they  can  be  commanded  over  inaccessible 
regions.  What  is  required  is  that  we  have  a sensor  located  at  a higher 
altitude  which  will  record  the  image  of  the  airborne  control  target  against 
the  natural  terrain  background  as  well  as  a Distance  Measurement  Navigation 
System  (DME)  to  continually  acquire  and  position  the  airborne  target 
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aircraft.  Both  sensor  and  navigation  systems  must  be  coupled  to  a 
precision  time  base. 

The  direct  ranging  or  distance  measurement  navigation  systems  operate 
in  the  microwave  frequency  region  1,000  to  15,000  Megacycles  which  limits 
observations  to  line  of  sight.  The  distance  measuring  equipment,  because 
of  short  wavelengths,  can  operate  in  either  a pulse  mode,  transmitting  a 
sequence  of  high  peak  power  sharp  rise  time  pulses  or  continuous  wave  mode, 
transmitting  three  or  more  frequencies  for  ambiguity  resolution.  Directional 
antennas  used  at  ground  locations  yield  high  si gnal -to-noi se  ratios  providing 
accuracies  in  the  order  of  a few  feet.  Directional  antennas,  because  of 
narrow  beams,  reduce  multipath  error  sources.  A pulse  DMF  system  is  pre- 
ferred for  the  application  in  this  report.  Large  scale  integrated  digital 
circuitry  produced  at  low  cost  provides  an  economical  means  of  implementing 
selective  coding  features;  and  low  cost  high  peak  power  transmitters  are 
available  for  airborne  transponders  making  the  pulsed  system  a preferred 
system  from  economical  and  flexibility  considerations.  The  high  accuracy, 
low  cost  and  compact  size  of  the  distance  measuring  system  has  opened  up  a 
number  of  applications  requiring  rapid  deployment  in  areas  of  sparse 
geodetic  control.  A common  requirement  of  most  applications  of  pulsed 
distance  measuring  equipment  is  the  accurate  location  of  multiple  ground 
positions  in  a network  tied  to  a geodetic  coordinate  system. 

An  experimental  high  precision  Distance  Measurement  System  was  recently 
installed  at  Rome  Air  Development  Center. (2)  The  description  of  this 
distance  measurement  type  is  included  as  Appendix  A,  "Distance  Measurement 
Equipment  Navigation". 


The  high  altitude  sensor  collects  aerial  photos  which  are  analytically 
reduced  based  on  the  general  bundle  adjustment  method.  An  aerial  photograph 
is  produced  by  a bundle  of  light  rays  that  pass  from  object  space  through 
the  camera  lens  to  the  film  in  the  image  plane.  At  the  instant  of  exposure 
image,  lens  and  object  occupy  discrete  points.  Knowledge  of  object  and 
image  positions  for  a finite  number  of  rays  enable  the  determination  of  the 
position  and  orientation  for  the  aerial  camera.  Utilizing  this  camera 
position  and  orientation  information,  one  is  able  to  measure  additional 
image  positions  and  compute  the  location  of  these  points  in  object  space. 

The  general  bundle  adjustment  method  refers  to  the  general  case  where 
the  number  of  photographs  and  the  number  of  rays  per  photo  may  vary  as 
circumstance  allows.  The  solution  is  a simultaneous  adjustment  of  all 
rays  on  all  photographs  to  produce  object  positions  for  all  points  of 
interest. 

Thus  the  DME  Navigation  System  provides  the  necessary  object  space 
control  to  enable  the  photo  collection  system  and  general  bundle  adjustment 
method  to  produce  space  positions  of  objects  not  under  the  control  of  the 
DME  Navigation  System. 

Appendx  B,  "General  Photogrammetric  Bundle  Adjustment  Solution", 
describes  the  mathematics  of  the  general  photogrammetric  bundle  adjustment 
solution. 

1.2  GENERAL  SYSTEM  CONCEPT  AND  PURPOSE  OF  THE  INVESTIGATION 

The  description  of  the  distance  measuring  network  and  remote  sensing 
aircraft  is  made  with  reference  to  Figure  1 and  Figure  2.  Figure  1 is  a 
typical  1:500,000  scale  map  showing  initially  installed  DME  ground 
navigation  stations  designated  20,  40,  70,  90  and  three  overlapping 
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photo  areas  tied  with  pass  points  10  through  15.  Aircraft  positions  are 
indicated  as  + points  1 through  9.  Shown  in  Figure  2 are  three  aircrafts 
positioned  by  DME  navigation  and  a remote  sensing  aircraft.  The  positions 
of  the  distance  measurement  equipped  aircraft  are  precisely  determined  by 
measurements  to  the  ground  stations  as  shown  in  Figure  1 which  anchors  them 
to  a geodetic  framework.  A remote  sensing  aircraft  or  satellite  with  a 
metric  camera  is  positioned  so  imagery  includes  the  DME  navigated  aircraft 
and  the  terrain  with  new  unknown  geodetic  positions  to  be  determined.  The 
photography  would  extend  into  a region  of  no  known  ground  control.  Employing 
analytic  aerotri angul ation  and  certain  known  parameters,  the  image  measure- 
ments of  unknown  locations  can  be  converted  to  known  ground  coordinates  in 
the  DME  coordinate  system.  Since  the  DME  position  is  electronically 
determined  and  tied  to  geodetic  coordinates,  all  features  in  the  imagery 
that  are  measured  with  respect  to  DME  aircraft  position  are  electronically 
determined  and  tied  to  geodetic  coordinates.  The  interface  between  the 
camera  equipped  remote  sensor  aircraft,  the  aircraft  under  DME  navigation 
control  and  the  DME  network  is  a precise  time  reference,  either  LORAN  C or 
the  clocking  system  in  the  DME  system  itself.  The  correlation  of  aerial 
photography  observations  with  DME  observations  on  the  same  imagery  enables 
all  the  observed  quantities  to  be  brought  to  a common  geometrical  base.  The 
solution  of  this  geometrical  problem  is  challenging  in  that  two  statistically 
independent  three  dimensional  systems  with  widely  varying  origin  locations 
and  orientations  are  superimposed  on  two  dimensional  photographs.  With  the 
aid  of  a precision  time  base  the  three  dimensional  relationships  are 


reconstructed. 


The  overall  purpose  of  the  study  is  to  set  up  a deterministic  model 
of  the  system  described  which  will  ^monstrate  its  feasibility.  Section 
2 describes  the  investigation  of  combined  DME  and  aerial  photography 
for  ground  control  extension.  Section  3 discusses  the  results,  and 
conclusions  and  reconmendations  are  given  in  Section  4. 
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2.0  INVESTIGATION  OF  COMBINED  DME  AND  AERIAL  PHOTOGRAPHY  FOR  GROUND 


1 


CONTROL  EXTENSION 

2.1  GENERAL 

The  object  of  the  main  experiment  in  this  study  is  to  test  the  concept 
of  extending  ground  control  with  the  aid  of  moving  airborne  targets.  The 
experiment  is  comprised  of  four  main  elements.  These  are  the  control  range, 
the  narrow  bandwidth  real-time  distance  measurements  from  the  control 
stations  to  the  moving  target  aircraft,  the  wideband  high  resolution  photo- 
graphy from  the  high  altitude  aircraft  which  images  the  moving  target 
aircraft  against  the  terrain  and  the  analytic  aerotriangulation  by  which 
the  ground  control  is  extended. 

A full  blown  test  of  the  procedure  was  not  deemed  practical  until 
preliminary  computations  demonstrated  the  feasibility  of  the  concept.  In 
order  to  make  the  computations  as  realistic  as  possible  an  existing 
instrumentation  range  at  Griffiss  AFB,  N.  Y.  was  utilized  and  the  experi- 
mental model  was  constructed  to  conform  to  the  dimensions  of  the  range. 

Bias  errors  and  noise  errors  are  measured  values  of  the  experiment. 

Section  2.2  describes  the  instrumentation  range,  Section  2.3  discusses 
the  distance  measurement  procedures  and  resulting  errors,  Section  2.4  covers 
the  formation  of  fictitious  photographs  to  adhere  to  the  instrumentation 
range  utilized  in  the  study,  and  Section  2.5  provides  details  on  the 
analytic  aerotriangulation  which  was  carried  out  to  extend  the  ground 
control  with  moving  airborne  targets.  The  results  are  listed  in  Section  2.6. 

2.2  INSTRUMENTATION  RANGE  AND  SURVEY  CONTROL 

The  reference  range  for  the  study  was  the  DME  Navigation  Instrumentation 
Sites,  RADC,  Griffiss  AFB,  New  York.  The  range  has  been  recently  surveyed 
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by  the  Geodetic  Survey  Squadron,  F.  F.  Warren  AFB,  Wyoming. ^ The  DMt 
Navigation  Range  is  illustrated  in  figure  1.  Four  control  stations  - 20, 

40,  70  and  q0,  were  used  in  this  study  and  only  they  are  shown  in  f igure  1. 
The  aircraft  positions  1 through  9,  as  shown  in  Figure  I,  were  presumed  to  be 
located  with  Distance  Measuring  fguipment  (DM1).  Ihe  master  stations  were 
located  at  stations  20,  40,  70  and  90. 

Table  1 lists  the  local  rectangular  space  coordinates  ot  the  above 
four  stations  together  with  the  estimated  achieved  survey  accuracies. 


TABU  1. 


SflfCTFO  SI  AT  IONS  IN  TUI  RADC  t'Ml  NAVIGATION 
RANG!  A1  GRIFF 1SS  Af  B , Nf W YORK 


local  Rectangular  Space  Coordinates  in  Meters 


Station 

X 

Y 

n 

l 

B 

20  (Tanner-ELS) 

527028.823 

♦ 0 . 3 

492266.378 

+0.2 

248.472 

+ 0 . 6 

40  (EP- 1 0337) 

1 

501204.070 

♦,0.7 

568050.961 

+0.3 

-24t> . 556 

+ 0.6 

70  (WP- 10339) 

479771 .957 

♦ 0.7 

56693 1 . 4t>5 

♦ 0.3 

- 72.468 

+0.6 

90  (RADC-ELS) 

479499. 707 

+0.1 

477793.045 

+0.3 

155.509 

+0.4 

2.3  DISTANCE  Ml ASURi Ml  NT  AND  SURVEY  fRRORS 

The  moving  airborne  target  procedure  for  control  extension  depends  on 
the  location  of  the  target  aircraft  at  the  instant  of  the  exposure  of  the 
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aerial  photograph  from  the  high  altitude  aircraft.  All  target  aircraft  and 
the  photography  aircraft  are  operating  on  a coimon  time  base.  The  error  in 
locating  the  target  aircraft  depends  upon  its  position  in  the  range  with 
respect  to  the  tracking  stations. 

Assume  that  survey  control  station  I is  tracking  aircraft  J.  At  a 
given  moment  when  the  high  altitude  aircraft  takes  a photo,  the  distance 
between  I and  J ( L -j  j ) is  : 

Li  j = C(Xj-Xi)2  ^ (Yj-Yi)2  + (Zj-Zi)23  ^ 0) 

Thus,  the  coordinates  of  the  target  aircraft  are: 

‘j  ■ CLU2-  (Vj-v,)2-  (zrzt )2J  s ♦ X( 

Vj  ■ Oij2  - (Xj-x,)2  - Uj-Z-i)2:  * V, 

Zj  « [Li/  - (Xj-Xi)2  - (Yj-Yi)2.]  ♦ Zi  (2) 

Errors  in  X j , Yj  and  Zj  are  due  to  the  errors  in  the  measured  distance 
L i i and  the  given  values  for  the  known  station  I ( X -j , Y i , Zi).  These  errors 
will  be  designated  ^ L i j , AXi , A Y i and  AZ,-. 

The  error  in  Xj  can  be  expressed  as  follows: 

ax  • ax.  ax.  ax. 

AXj  = — J-  Alii  + — 1 AXi  + — J-  AY i + — 1 AZi  (3) 

3Lij  3Xi  aY i aZi 
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Under  the  assumption  of  no  correlation  between  the  variables  and  the 
equivalence  of  °Xj,  °X^,  °Y^  and  to  aXj,  aX.,- , aY^  and  AZ^ 
respectively: 


°X; 


3XH 


3L  • • 
TJ 


o 3Xi  p 9 X n 3 X j 

OL,,:2  + O-^oXiT  + [—  aYi]2  + c-2  02,] 
'J  3X-  av.  S7. 


(4) 


3Y; 


9Zi 


In  a similar  fashion  expressions  for  °Yj  and  °Zj  can  be  generated.  °X.j , 
oy..  and  oZ-  values  were  taken  from  Table  1 which  are  the  results  of  an 
actual  survey.  DME  measurements  at  the  RADC  Navigation  Range  place  the  oL^j 
value  at  3.14m.  The  DME  navigation  system  is  described  in  Appendix  A. 

Utilizing  the  RADC  DME  Navigation  Range  and  stations  20,  40,  70  and  90 
as  shown  in  Figure  3 an  error  propagation  study  was  conducted  for  the 
target  aircraft  located  at  each  of  the'  range  positions  (+).  The  value 
6m.  at  position  Xj  = 53,  Yj  = 53  was  the  average  of  four  Xj  values 
computed  from  stations  20,  40,  70  and  90. 

Figure  4 expresses  the  reliability  of  the  values  given  in  Figure  3. 
Similarly,  Figures  5 and  6 express  the  error  propagation  in  Yj  and  Figures 
7 and  8 apply  to  the  Zj  coordinate  of  the  target  aircraft. 

The  resultant  errors  are  displayed  in  Figure  9 along  with  standard 
deviation  values. 

where: 

Resultant  = [Xj2  + Yj2  + Zj^]15  (5) 
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It  should  be  noted  that  70  to  80  percent  of  the  error  values  in  Figure 
9 are  due  to  the  error  propagation  in  Z.  The  planimetric  errors  over  this 
area  of  the  range  would  average  about  10  meters. 

2.4  HIGH  ALTITUDE  PHOTOGRAPHY 

The  general  concept  of  the  use  of  high  altitude  photography  of  moving 
airborne  targets  against  the  natural  terrain  background  is  shown  in  Figure 
2.  The  high  altitude  aircraft  exposes  a photo  and  images  two  types  of 
points.  The  first  is  the  moving  airborne  targets  for  which  X,  Y,  and  Z 
are  known.  The  second  is  pass  points  whose  X,  Y,  and  Z coordinates  must 
be  computed. 

In  this  study  three  fictitious  high  altitude  photographs  were  created. 
Figure  1 illustrates  the  position  and  extent  of  coverage  of  these  photos 
within  the  RADC  DME  Navigation  Range. 

Tables  2,  3 and  4 summarize  the  parameters  of  the  three  fictitious 
high  altitude  photographs. 

2 . 5 ANALYTIC  AEROTRI ANGULAT ION  COMPUTATIONS 

Utilizing  the  RADC  DME  Navigation  Range  and  a standard  bundle  adjust- 
ment solution,  a fictitious  flight  line  of  three  consecutive  aerial 
photographs  was  triangulated  for  the  solution  of  six  unknown  ground  points 
and  the  three  exposure  stations.  Figure  1 shows  the  location  of  the  photos 
within  the  range.  Recall  that  points  1 through  9 are  aircraft  which  are 
located  by  DME  and  appear  on  the  photos  as  targets.  Points  10  through  15 
are  selected  pass  point  images  whose  range  coordinates  will  be  determined 
as  a result  of  the  analytic  aerotriangulation. 

The  input  to  the  solution  includes  the  DME  determined  aircraft 
coordinates  and  the  photo  coordinates  of  the  aircraft  and  selected 
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TABLE  2.  FICTITIOUS  HIGH  ALTITUDE  PHOTO  NO.  1 


X]c=500,500  m. 

Xp]=0  mm. 

Y]c=527,400  m. 

ypl=0  irm. 

Z]c=20,000  m. 

f=152.4  run. 

K=0°30 1 

+ - Moving  Airborne  Targets 

}s  l°30 ' 

0 - Pass  Points 

W=0°30 1 

f - Focal  Length 

*1 

-78.185* 

*2 

74.977 

x3  3.522  x 1 g 

+3.000 

VI 

98.988 

yz 

100.790 

y3  -53.749  yiO 

+100.000 

X1 

495,000.000**  X2 

505,000. 

000  X3  500,500.000  X10 

500,256.576 

534,000.000  Y2 

534,000. 

000  Y3  524,000.000  Y10 

540,643.380 

Z1 

10100.000 

z2 

10150.000  Z3  10,075.000  Z10 

200.000 

xn 

1 04 . 000 

x12 

2.000  x-j  3 98.000 

yn 

105.000 

*12 

-90.000  y 1 3 -95.000 

xn 

513,322.474 

x12 

500,344.655  X13  512,637. 

071 

vn 

541 ,283.620 

Y12 

515,972.132  Y13  515,521.093 

Zll 

50.000 

z12 

250.000  Zl3  60.000 

* 

- Photo 

Coordinate  in  mm. 

- Range  Coordinate  in  m. 
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TABLE  3.  FICTITIOUS  HIGH  ALTITUDE  PHOTO  NO.  2 


+y 


+x 


X2C=512,000  m. 
y2c=528,500  m. 
Z2c=20 ,000  m. 
K=0°30 ' 

»=1° 

W=0°45 ' 


Xp2=0.0  mm. 

yp2=0.0  mm. 
f = 1 52 .4  mm. 

+ - Moving  Airborne  Targets 
0 - Pass  Points 
f - Focal  Length 


x4 

-78.189* 

x5 

96.194 

x6 

9.572 

*4 

87.734 

V5 

100.603 

X6 

-86.401 

x4 

506,500.000** 

x5 

518,000.000 

x6 

512,500.000 

n 

534,500.000 

y5 

535,200.000 

Y6 

523,000.000 

u 

10,100.000 

z5 

10,310.000 

z6 

10,175.000 

x10 

-84.480 

xn 

13.415 

x12 

-87.083 

yio 

89.674 

yn 

93.919 

yi2 

-97.762 

x10 

500,256.576 

xn 

513,322.474 

X1 2 

500,344.655 

Y10 

540,643.380 

Yn 

541 ,283.620 

Y12 

515,972.132 

z10 

200.000 

zn 

50.000 

Zl2 

250.000 

X1 3 

6.644 

x14 

99.000 

X1 5 

90.000 

yi3 

-101.193 

yu 

93.000 

yi5 

-85.000 

X1 3 

512,637.071 

X14 

524,372.360 

X1 5 

523,302.770 

Y1 3 

515.521.093 

y14 

540,918.912 

Y1 5 

517,951.376 

Z1 3 

60.000 

z14 

400.000 

z15 

300.000 

★ - 

★ ★ _ 

Photo  Coordinate  in  mm 
Range  Coordinate  in  m. 
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TABLE  4.  FICTITIOUS  HIGH  ALTITUDE  PHOTO  NO.  3 

X3C=523,000  m.  Xp3=0.0  mm. 

Y3C=529,000  m.  yp3=0.0  mm. 

Z3C=20,100  m.  f=152.4  nm. 

K=l°  + - Moving  Airborne  Targets 

<p=2°  0 - Pass  Points 

W=l°  f - Focal  Length 


*7 

-85.000 

x8 

+92.000 

x9 

+2.000 

xn 

-64.811 

*7 

+88.000 

n 

+90.000 

y9 

-101 .000 

yn 

89.474 

x7 

516,787,773  X8 

528,472.625 

x9 

522,899,415 

xn 

513,322.474 

Y7 

535,029.936  Y8 

535,067.862 

* Y9 

522,657.776 

Yn 

541 ,283.620 

ll 

10,100 

z8 

10,200 

z9 

10,150 

Zll 

50.000 

X1 3 

-74.912 

x14 

17. 

416 

X1 5 

6.131 

Y13 

-103.276 

yi4 

88. 

587 

yi5 

-88.786 

X13 

512,637.071  X14 

524 

,372. 

,360  X15 

523,302, 

.770 

Y1 3 

515,521.093  Y14 

540 

,918.912  Y15 

517,951 , 

.376 

Z13 

60 . 000 

Z14 

400.000 

Z15 

300.000 

* - Photo  Coordinate  in  mn. 

**  - Range  Coordinate  in  m. 
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pass  points.  Since  these  are  measured  quantities  a weight  is  assigned  to 
each  of  the  values.  Approximations  are  assigned  to  the  unknown  quantities, 
namely,  the  outer  orientation  values  and  the  unknown  range  coordinates  for 
the  pass  points.  Weights  are  used  as  constraints  in  the  solution  and 
accordingly  the  approximate  values  are  given  a weight  of  zero  so  that  they 
may  have  complete  freedom  in  the  solution  and  obtain  their  most  probable 
values . 

Twelve  aerotriangulation  solutions  were  computed.  The  following  four 
had  particular  significance. 

1)  Solution  I - The  purpose  of  this  solution  was  to  demonstrate  that 
the  deterministic  aerotriangulation  model  and  computational  procedure  were 
correct.  All  input  control  and  coordinate  information  are  correct  and  the 
solution  converges  to  the  correct  values  for  the  outer  orientation  para- 
meters and  the  range  coordinates  for  the  pass  points. 

2)  Solution  2 - In  this  solution  the  range  coordinates  of  the  target 
aircraft  were  perturbed  according  to  the  error  propagation  resulting  from 
DME  and  ground  station  survey  errors.  The  appropriate  X,  V and  Z errors 
were  selected  from  Figures  3,  5 and  7.  Signs  for  these  error  values  were 
selected  from  a table  of  random  numbers.  A weight  of  unity  was  used  for 
the  above  coordinates.  The  photo  coordinates  were  given  the  weight  of 
2500  in  this  solution.  This  assumes  a 20  micrometer  standard  deviation  in 
each  photo  coordinate. 

3)  Solution  3 - In  this  solution  the  range  coordinates  of  the  target 
aircraft  were  perturbed  in  the  same  manner  as  in  solution  2.  However,  in 
this  solution  the  weights  used  for  the  above  coordinates  reflect  the  error 
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introduced  into  the  coordinate.  The  weight  of  the  X ^ordinate  of  target 
aircraft  1 is  0.02  since, 

W,  = -L-  = _L  = .02 

°X1 2 72 

thus  reflecting  the  7 meter  error  introduced  into  X] . 

The  photo  coordinates  were  given  the  weight  of  2500  in  this 
solution.  This  assumes  a 20  micrometer  standard  deviation  in  each  photo 
coordinate. 

4)  Solution  4 - This  solution  was  similar  to  3 except  that  the  photo 
coordinates  in  this  case  were  weighted  500  rather  than  2500.  Thus,  it  was 
assumed  that  each  photo  coordinate  may  have  a standard  deviation  of  45 
mi crometers . 

2.6  RESULTS  OF  ANALYTIC  AEROTRIANGULATION 

The  purpose  of  this  study  was  to  demonstrate  how  moving  airborne 
targets  could  be  utilized  to  extend  ground  control.  The  results  are 
presented  in  the  form  of  a number  of  figures  showing  the  area  of  the  range 
that  control  was  being  extended  into  and  the  actual  errors  that  resulted 
from  the  aerotriangulation  as  previously  described.  Error  is  defined  as 
the  absolute  value  in  meters  between  the  known  coordinate  from  the  deter- 
ministic model  and  the  final  coordinate  values  given  by  the  analytic 
aerotriangulation  solution. 

The  results  of  solution  2 are  illustrated  in  Figure  10.  The  values 
in  the  parenthesis  are  the  absolute  resultant  errors  that  apply  to  the 
nearest  pass  point  or  exposure  station.  Resultant  is  defined  in  equation 
(5).  Xr  equals  the  average  resultant  error  to  the  nearest  meter  for  the 
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nine  points,  namely;  the  six  pass  points  and  the  three  exposure  stations. 
or  equals  the  standard  deviation  to  the  nearest  meter  for  these  nine  values. 
Contours  have  been  interpolated  based  on  the  nine  discrete  point  values. 

Values  for  solution  3 are  shown  in  Figure  11.  Once  again  the  errors 
are  resultant  values. 

The  most  favorable  results  were  obtained  from  solution  4.  Figures  12, 
13,  14,  and  15  are  used  to  display  the  results  of  solution  4.  Figure  12 
gives  the  resultant  errors  and  Figures  13,  14,  and  15  show  the  X,  Y,  and 
Z errors  respecti vely. 
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3.0  DISCUSSION  OF  RESULTS 


The  results  will  be  discussed  in  terms  of  three  questions.  First, 
what  errors  could  one  expect  in  the  location  of  ground  control  points  in 
the  range  if  DME  and  control  survey  errors  were  the  error  sources;  second, 
what  errors  could  one  expect  in  the  location  of  ground  control  points  in 
the  range  if  analytic  aerotriangulation  were  utilized;  and  lastly,  how  do 
the  two  techniques  compare? 

Using  actual  error  values  from  DME  experience  and  surveys  conducted 
on  the  RADC  DME  Navigation  Range,  Figure  16  summarizes  the  error  propagation 
values  one  would  expect  for  the  various  locations  within  the  range. 
Considering  the  49  discrete  points  sampled  the  average  resultant  error 
to  the  nearest  meter  is  25m.  and  the  standard  deviation  to  the  nearest 
meter  is  2m.  As  is  shown  in  Figure  7,  the  largest  component  of  the  25m. 
error  is  the  Z value  which  accounts  for  approximately  20m.  of  the  error. 

Figures  10,  11  and  12  summarize  the  results  of  the  three  solutions  of 
interest  for  the  analytic  aerotriangulation.  The  results  in  Figure  11 
(solution  3)  are  superior  to  those  of  Figure  10  (solution  2).  The  values 
in  solution  3 are  more  consistent  as  the  standard  deviation  of  a single 
point  is  reduced  from  11m.  to  8m.  Thus  the  weighting  of  the  control  to 
real istical ly  reflect  the  errors  introduced  yields  a better  solution. 

Continuing,  it  is  seen  that  Figure  12  (solution  4)  improves  over 
Figure  11  (solution  3).  Solution  4 utilized  a photo  coordinate  weighting 
which  allows  for  an  error  of  45  micrometers  in  the  photo  coordinates  rather 
than  the  20  micrometers  that  was  allowed  in  solution  3.  in  effect  the 
photo  coordinates  in  4 were  allowed  more  flexibility  in  seeking  the  best- 
fit  solution.  Under  the  dynamic  conditions  for  the  moving  airborne 
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targets,  the  45  micrometer  flexibility  in  photo  coordinates  may  very  well 
be  more  realistic  than  20  micrometers. 

The  improvement  in  solution  4 is  both  in  average  resultant  error  and 
the  standard  deviation  value. 

A comparison  of  Figures  12  and  16  will  give  one  an  idea  of  the  relative 
merits  of  analytic  aerotriangulation  control  extension  versus  the  use  of  the 
DME  solutions.  Figure  17  summarizes  this  comparison.  Based  on  this  study 
the  analytic  aerotriangulation  is  superior  to  DME  alone  in  the  shaded  areas 
of  Figure  17. 

It  is  of  interest  to  note  that  the  largest  error  component  for  the 
analytic  aerotriangulation  solution  is  the  Y coordinate.  The  errors  in  X 
are  less  since  the  control  for  the  solution  was  geometrically  stronger  in 
the  X direction.  The  reverse  was  true  for  the  DME  solutions  as  was  shown 
previously.  The  aerotriangulation  solutions  for  Z were  greatly  improved 
over  DME  solutions  as  one  would  expect  under  the  greatly  improved  base/ 
height  ratios  existing  in  the  aerotriangulation  case. 

There  appears  to  be  no  systematic  reasons  for  the  differences  as  shown 
in  Figure  17.  Keep  in  mind  that  all  error  values  and  sums  are  absolute 
without  defined  direction. 

It  is  interesting  to  consider  the  moving  airborne  target  concept  from 
a practical  viewpoint.  The  study  assumes  that  four  aircraft  were  in  the 
air  simultaneously  and  that  there  were  a minimum  of  four  DME  ground  stations 
acquiring  and  locating  the  airborne  control  targets  on  a precisely  coordinated 
time  base.  This  is  admittedly  an  elaborate  system,  however,  it  does  result 
in  the  location  of  ground  stations  without  occupying  the  stations  or,  in 
fact,  being  able  to  directly  view  them.  DME  solutions  alone  would  require 
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that  the  stations  be  occupied.  In  addition,  all  currently  occupied  stations 
and  stations  to  be  occupied  must  be  intervisible.  The  moving  airborne 
targets  have  considerable  advantage  here  since  1 ine-of-sight  problems  are 
less  severe. 

The  system  can  be  simplified  by  using  only  one  moving  airborne  target 
which  flies  systematical ly  about  the  area  of  interest.  This  single  air- 
craft would  be  photographed  by  the  high  altitude  aircraft  numerous  times 
and  again  the  analytic  aerotriangulation  solution  would  be  employed. 

The  timing  and  DME  aerotriangulation  correK  ‘•ion  problems  seem  complex 
but  the  state-of-the-art  is  moving  very  rapidly  in  this  area  and  it  seems 
that  systems  of  this  nature  will  be  routinely  handled  in  the  future. 

The  old  problems  of  image  measurement  and  point  identification  are 
still  with  us  in  this  application.  There  doesn't  seem  to  be  any  reason 
for  them  to  be  more  severe  in  this  case,  however,  only  an  actual  test  of 
the  procedure  will  answer  this  question. 


4.0  CONCLUSIONS  AND  RECOMMENDATIONS 


The  following  conclusions  have  been  drawn  from  the  study. 

1)  Narrowband  distance  measurements  have  been  successfully  combined 
with  wideband  photography  for  the  extension  of  range  control.  There  is  an 
indication  that  the  concept  may  be  as  good  as  standard  ground  techniques 
for  ground  control  extension. 

2)  The  moving  airborne  target  technique  is  very  flexible  and  photo 
flight  lines  can  be  run  which,  due  to  superior  geometry  and  redundant 
sampling,  provide  small  errors  in  X,  Y and  Z coordinates.  The  flight  line 
orientation  in  the  subject  study  favored  the  X coordinate.  Command  of  the 
aircraft  flight  path  to  selected  areas  can  minimize  errors.  The  desired 
geographical  flight  areas  can  be  determined  prior  to  aircraft  flight. 

3)  The  procedures  described  represent  an  order  of  magnitude  decrease 
in  the  time  required  to  perform  ground  surveys  to  DME  navigation  system 
requi remen ts . 

4)  The  feasibility  of  the  procedure  will  improve  as  each  of  the 
system  components  approaches  a real  time  capability.  This  includes  real- 
time data  reduction  by  the  DME,  together  with  the  use  of  high  resolution 
charge  coupled  device  cameras  programmed  for  real  time  solutions. 

The  following  recommendations  are  made: 

1)  Additional  effort  needs  to  be  expended  to  determine  the  most 
efficient  use  of  the  moving  airborne  target  concept. 

2)  This  technique  applies  a combination  of  DME  and  photo  sensing 
instruments.  Similar  concepts  using  a combination  of  other  sensors  should 
be  explored.  In  this  way  each  sensor  can  contribute  its  outstanding 
characteristic  to  the  total  mission. 
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APPENDIX  A 


DISTANCE  MEASUREMENT  EQUIPMENT  (DME)  NAVIGATION 

This  appendix  describes  one  kind  of  ground-air  Distance  Measurement 
Equipment  (DME)  Navigation  System,  using  pulsed  data,  that  meets  the 
precision  required  for  extension  of  geodetic  position.  In  addition  to 
determining  accurate  position  and  derivatives  of  position,  the  link 
transmits  data,  includes  error  correcting  codes  (Hanning-Reed  Solomon), 
and  is  a full  duplex  system  operating  at  L band. 

Figure  A-l  is  a simplified  diagram  showing  the  basic  airborne 
portion  of  the  Distance  Measurement  Equipment  and  the  ground  station. 

A master  navigation  computer  is  located  at  a Ground  Control  Station  and 
contains  a Kalman  filter  implementation  that  performs  A/C  to  grid, 
coordinate  transfer,  performs  integration  of  measurements,  does  pre- 
filtering of  airborne  measurements  (when  air  derived  sensors  are  used) 
and  bias  error  correction.  The  transmitted  messages  that  includes  a 
ranging  pulse  are  made  up  of  8 words  extending  over  500  usees.  Each 
word  consists  of  a preamble,  eight  triple  pulses  and  a ranging  pulse  at 
the  end  of  the  sequence.  Figure  A-2  illustrates  the  format  of  transmitted 
information.  The  initiation  of  the  Radio  Frequency  (RF)  transmission  sets 
a clock  at  the  DME  Ground  Control  Station  (GCS)  and  the  reception  of  the 
signal  at  aircraft  and  from  the  DME  GCS  through  the  aircraft  to  ground 
beacon  stations  sets  clocks  at  these  locations.  The  range  is  determined 
by  comparing  clock  times.  These  clock  times  obtained  from  various  trans- 
mission paths  are  converted  to  range  values  in  the  DME  GCS  navigation 
computer. 


A-l 


Simplified  transmission  paths  for  a GCS  and  two  remote  beacon/ 
transponder  configurations  are  shown  in  Figure  A-3.  Path  1 is  from  the 
GCS  to  Aircraft  No.  1,  through  the  airborne  hardware,  calibrates  the 
equipment  delays,  initiates  the  aircraft  clock  and  retransmits  to  the  GCS. 
Path  2 goes  to  the  Aircraft  No.  1 and  is  redirected  to  Remote  Beacon  1 
where  the  beacon  hardware  delays  are  calibrated,  the  beacon  clock  is  reset 
and  the  message  is  retransmitted  to  the  GCS  through  the  Aircraft  No.  1. 

The  range  measurement  is,  therefore,  the  round  trip  path  from  the  Beacon 
No.  1 to  the  Aircraft  No.  1 and  hence  to  the  GCS.  The  sequence  to  obtain 
calibration  of  hardware  delays,  send  messages  to  three  aircraft  and  two 
ground  beacons,  position  location  requires  transmission  times  on  the  order 
of  4 milliseconds. 

The  random  error  in  range  is  determined  by  referencing  time  to  the 
half  power  point  on  the  leading  edge  of  the  last  pulse  in  the  message 
sequence  in  Figure  A-2.  The  variance  in  a single  observation  is  directly 
related  to  the  slope  of  the  leading  edge  and  inversely  proportional  to  the 
square  root  of  the  signal  energy  to  noise  power  and  the  number  of  single 
observations  processed  prior  to  a range  determination.  Single  range 
observations  are  made  at  the  rate  of  30  per  second.  Ten  measurements  are 
integrated  prior  to  a range  readout.  Random  errors  in  time  determination 
on  a single  observation  are  on  the  order  of  10-15  nanoseconds  for  20  dB 
signal-to-noise  ratio.  Residual  bias  errors  are  7-15  nanoseconds.  Certain 
error  sources  affecting  the  random  and  systematic  errors  have  corrections 
applied. 

Multiple  path  effects  are  eliminated  by  correlating  Triplets 
received  with  known  message  structure.  Standard  corrections  are  made 
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for  propagation.  Angles  of  transmission  are  normally  high  virtually 
illuminating  residual  propagation  errors,  and  the  frequency  of  operation 
is  optimum  with  regard  to  environmental  noise.  Major  component  of  noise 
is  the  electronic  equipment  and  variation  in  signal  energy  due  to  range 
dependent  (R?)  loss  and  nulls  in  antenna  patterns. 

A simplified  diagram  of  the  computational  process  is  shown  in  Figure 
A-4.  Computation  is  done  at  the  GCS.  Computation  includes  the  co-variance 
matrices  and  Kalman  filter  states  for  each  aircraft.  Serial  processing  is 
performed  on  the  navigation  data.  The  navigation  data  includes  three 
position,  velocity  and  acceleration  for  each  aircraft  referenced  to  the 
geocentric  inertial  frame.  Computation  is  done  in  the  aircraft  frame 
of  reference  (DME  information  is  also  transformed  into  the  platform  frame). 
The  Kalman  filter  has  twelve  states  as  initial  conditions;  they  are  3 
axis  inertial  platform  misalignment,  3 aircraft  position  states,  3 aircraft 
velocity  states  and  3 inertial  platform  gyro  states.  Only  the  x,  y,  z, 
x,  y,  z for  three  aircraft  are  of  concern;  therefore,  the  navigation  and 
measurement  error  co-variance  matrices  are  18  x 18.  Data  on  the 
effectiveness  in  reducing  range  and  velocity  errors  through  complimentary 
filtering  has  yet  to  be  obtained. 

Emphasis  is  on  improving  aircraft  velocity  rather  than  position. 
Therefore,  position  measurement  originally  submitted  should  be  used  until 
measurements  confirm  a change  is  warranted.  Of  course,  Kalman  filtering 
cannot  reduce  bias  error. 
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GENERAL  PHOTOGRAMME TR I C BUNDLE  ADJUSTMENT  SOLUTION 

1 . Introduction 

The  procedure  to  be  employed  is  a simultaneous  photogrammetric  bundle 
adjustment  solution.  All  computations  concerning  photo  (sensor)  orientation 
and  position,  DME  control  coordinates,  photo  coordinates  and  added  para- 
meters are  handled  with  a unified  approach.  It  is  defined  as  the  least 
square  solution  for  all  unknown  parameters  using  all  of  the  image  and  DME 
Navigation  Range  coordinate  observations  for  the  photos  in  question. 

This  procedure  is  highly  flexible  as  it  can  be  applied  to  photographs, 
digital  records,  blocks  of  photography,  etc.  The  actual  degree  of 
flexibility  will  obviously  be  constrained  by  the  cost-effectiveness  of 
the  computations. 

The  following  discussion  briefly  describes  the  general  photogrammetric 
bundle  adjustment  solution  (simultaneous  record  adjustment).  Paragranh  2 
covers  the  basic  condition  equation.  The  linearized  condition  is  given  in 
Paragraph  3.  General  and  special  solutions  are  discussed  in  Paragraph  4. 

2.  Basic  Condition  Equation 

The  collinearity  equation  is  the  basic  condition  equation  which  relates 
the  DME  range  coordinates  to  the  photo  coordinates.  This  equation  type  is 
the  basis  for  single  photo  resection,  relative  orientation  and  strip  or 
block  aerotriangulation. 
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coordinates  are  related  to  photo  coordinates  in  these  equations. 


xij  ' dxi j = _ci 

Vij  ‘ dyij  = 'ci  — 
A3 

where 


Al 

1 

U 

X 

r~i 

X 

1 

A 2 

= Mi 

Yj  - Yi c 

a3 

Zj  - ZiC 

Equation  terms  are 


(1) 


(2) 


Ci  - principal  distance  for  i-th  photo. 

xij,yij  - photo  coordinates  for  j-th  point  of  i-th  photo. 

dxij.dyij  - distortions  to  ij-th  image  because  of  any  added 
parameters 

M.j  - orthogonal  rotation  matrix  for  i-th  photo. 


Mi 

where  C0  = 


C0CK  : CWSK  + SWS0CK  : 

SWSK 

- SWS0CK 

-C0CK  : CWCK  - SWS0SK  : 

SWCK 

+ CWS0SK 

(3) 

S0  : -SWC0  : 

CWC0 

L 

** 

cos  phi;  SK  = sin  kappa;  CW 

= cos 

omega;  etc 

. Omega, 

phi  and  kappa  are  primary,  secondary  and  tertiary  rotations. 


Xj.Vj.Zj  - ground  coordinates  of  j-th  point. 

Xic,Y.jC,Zic  - ground  coordinates  of  camera  station  for  i-th  photo. 
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3.  Linearized  Conditions 

The  basic  collinearity  condition  is  nonlinear.  Since  most  solutions 
will  use  a least  squares  adjustment  and  these  computations  are  generally 
performed  with  linear  functions,  it  is  necessary  to  linearize  the  basic 
conditions.  Series  expansions,  and  Taylor  Series  in  particular,  are 
generally  used  for  this  process  where  only  the  zero  and  first-order  terms 
are  included. 

The  linearized  expression  for  the  collinearity  equations  can  be 
presented  in  matrix  notation  as  follows: 


oo  p p G G 

vij  + Bij  sij  + Bij  s + Bij  sj  + eij  = 0 

( a , 1 ) (a ,b) ( b, 1 ) (a,c)(c,l)  (a ,3) ( 3,1 ) (a,l) 


where 


* vector  of  residuals  for  the  observation  of  the  ij-th  image. 


matrix  of  partial  derivatives  for  the  orientation  parameters 
of  the  i-th  unit. 


matrix  of  linear  coefficients  for  the  added  parameters. 


mat'-'x  of  partial  derivatives  for  the  ground  coordinates  of 
the  j-th  point. 


vector  of  corrections  for  the  orientation  parameters  of  the 
i-th  unit. 


vector  of  added  parameters. 


vector  of  corrections  for  the  ground  coordinates  of  the 
j-th  point. 
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1 


= discrepancy  vector  for  the  ij-th  image. 

= number  of  observation  equations  per  photogrammetric  point 
that  is  observed.  (2  for  photos  - x,y  coordinates) 

= number  of  orientation  parameters  per  unit  (6  for  each  photo). 

= number  of  added  parameters. 


If  "a  priori"  values  are  available  for  the  orientation  parameters, 
added  parameters,  or  ground  coordinates  they  can  be  treated  as  observations 
as  follows: 


0 

vi 

- si 

+ ei 

= 0 (orientation  parameters  of 

the  i-th  unit) 

(5) 

(b,l)(b,l)(b,l) 

p 

p 

p 

V 

- s 

+ c 

= 0 (added  parameters) 

16) 

(c,l)(c,1)(c,l) 

G 

G 

G 

VJ 

- SJ 

= 0 (ground  coordinates  of  the 

j-th  point) 

(7) 

(3, 

1)(3,1)(3,1) 

where 

0 

p 

G 

vi’ 

V. 

vj 

= vectors  of  residuals  for  the 

observations . 

0 

p 

G 

ei » 

c , 

j 

= discrepancy  vectors  for  the  observations. 

The  weight  matrices  for  the  observations  are: 

w-jj  = image  observations 
(a, a) 


orientation  parameters 


o 

wi 

(b.b) 

p 

w = added  parameters 

(c,c) 

G 

wj  = ground  coordinates 

(3,3) 

4.  General  or  Special  Solutions 

The  completely  general  solution  would  allow  for  the  simultaneous 
solution  of  equation  types  (4),  (5),  (6)  and  (7).  Special  solutions  would 
allow  for  all  combinations  of  the  above  equation  types.  The  weight  matrices 
would  be  used  as  constraints  allowing  various  parameters  to  be  treated  as 
known  or  unknown. 

As  an  example,  suppose  a stereo  pair  were  being  evaluated  and  there 
were  ten  points  of  interest  in  the  stereo  model.  Of  these  ten  points,  the 
ground  coordinates  were  known  for  4 of  them.  The  orientation  parameters 
for  each  photo  are  unknown,  5 added  lens  distortion  parameters  are  unknown 
and  the  X,  Y,  Z coordinates  of  six  of  the  object  points  are  unknown.  For 
this  solution,  the  following  2 types  of  linearized  observation  equations 
are  to  be  used: 

oo  p p G G 

vij  + Bij  Si  + Bij  s + Bij  sj  + £ij  = 0 (g) 
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0 


(9) 


VJ  ■ Sj  + £J 


In  this  case  "a  priori"  values  are  provided  only  for  the  ground 
coordinates  of  the  four  known  points.  All  other  unknowns  will  be  approxi- 
mated and  weights  of  zero  supplied  in  each  case. 

In  this  example,  there  will  be  40  type  (8)  equations  and  30  type  (9) 
equations. 

These  70  equations  are  expressed  as  follows: 

v + BS  + r=  0 (1C 


where 
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The  reduced  normal  equations  will  be  of  rank  (6xm+3xn+c) 
where  m equals  the  number  of  photos,  n equals  the  number  of  ground  points 
and  c equals  the  number  of  added  parameters.  For  the  given  example  there 
will  be  [(6x2)  + (3x10)  + 5 = 47]  normal  equations,  one  for  each  unknown. 
The  70  observation  equations  are  now  reduced  to  the  47  normal  equations 
which  are  symbolized  as  follows: 


N S + C = 0 

(47,47) (47,1 ) (47,1) 


where 


0 


_T  _ 

N = B W B ; C 


B W e ; W 


P 

0 W 0 


G 

0 0 W 


(11) 


I 


The  solution  for  corrections  to  first  approximations  would  be  calculated 
as  follows: 

-1 

S = -N  C (12) 

The  ? quantities  would  be  added  to  the  first  approximations  and 

sufficient  iterations  processed  for  the  sum  of  the  weighted  residuals 

__T 

squared  to  become  a minimum  (V  W V -*  minimum).  Thus,  the  best  estimate1- 
for  the  unknowns,  47  in  this  case,  are  obtained. 


8-7 


as  a by  product 


The  precision  of  each  calculated  unknown  is  obtained 
of  the  least  squares  solution. 

Covariance  matrix  U)  = ^ n'1 

where 

'2  _J_  _ 

a0  = V w v 

y 

v - -7 

V = exampl'e°f  degreeS  °f  freedom  which  is  23  in  the  preceding 

The  standard  deviations  of  the  unknowns  are  the  square  roots  of  the 
diagonal  terms  of  the  covariance  matrix. 

The  general  solution  as  modified  by  weight  constraints  can  be  used 
in  many  combinations,  one  of  which  has  been  used  as  an  example. 
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MISSION 

of 

Rome  Air  Development  Center 

RAVC  plant  and  execute*  research,  development,  test  and 
& elected  acquisition  programs  in  support  of  Command,  Control 
Communications  and  Intelligence  (C3I)  activities.  Technical 
and  engineering  support  within  areas  of  technical  competence 
Is  provided  to  ESI)  Program  0$  faces  (PC*)  and  otheA  ESV 
elements.  The  principal  technical  mission  areas  axe 
communications,  electromagnetic  guidance  and  control,  sur- 
veillance of  ground  and  aerospace  objects,  intelligence  data 
collection  and  handling,  information  system  technology, 
ionospheric  propagation,  solid  state  sciences,  microwave 
physics  and  electronic  reliability,  maintainability  and 
compatibility. 


